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Abstract 

We construct a nearest-neighbor Hamiltonian whose ground states 
encode the solutions to the NP-complete problem INDEPENDENT 
SET in cubic planar graphs. The Hamiltonian can be easily simulated 
by Ising interactions between adjacent particles on a 2D rectangu- 
lar lattice. We describe the required pulse sequences. Our methods 
could help to implement adiabatic quantum computing by "physically 
reasonable" Hamiltonians like short-range interactions. 

1 Introduction 

Adiabatic quantum computation has been proposed as a general way 
of solving computationally hard problems on a quantum computer 
FeaOl . Adiabatic quantum algorithms proposed so far work by ap- 
plying a time-dependent Hamiltonian 

H(t) = (l-t)H B + l Hp (1) 

that interpolates linearly from an initial Hamiltonian Hp to the final 
Hamiltonian Hp. The Hamiltonians are chosen such that the ground 
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states of Hp are easily prepared and the ground states of the final 
Hamiltonian Hp encode the solutions to the problem |Fea01| . 

The running time of the algorithm is denoted by T. If H (t) varies 
sufficiently slowly, i.e., T is sufficiently high, then one hopes that the 
final state of the quantum computer will be close to the ground state 
of the final Hamiltonian Hp, so a measurement will yield a solution 
to the problem with high probability. The adiabatic theorem is the 
justification for this hope. However, it is not clear whether all neces- 
sary conditions for adiabatic evolution are satisfied. For instance, it is 
not clear whether the gap between the ground states and first excited 
states of H(t) is sufficiently high for all t. 

The adiabatic method can only succeed if the Hamiltonian H{t) 
changes slowly. But how slow is slow enough? Unfortunately, this 
question has proved difficult to analyze in general. Some numerical 
evidence suggests the possibility that the adiabatic method might ef- 
ficiently solve computationally interesting instances of hard combina- 
torial problems, outperforming classical algorithms FeaOl . Whether 
adiabatic quantum computing provides a definite speedup over classi- 
cal methods for certain problems remains an interesting open question. 

Our objective in this paper is not to explore the computational 
power of the adiabatic quantum computing, but rather to investigate 
how to implement the adiabatic time evolution starting from "physi- 
cally reasonable" Hamiltonians like short-range interactions. 

A Hamiltonian can be considered as physically reasonable only if 
it is "local". One way to describe locality is a follows. Let TL := C 2 
denote the Hilbert space of a single qubit and the joint Hilbert 
space of n qubits. L(W® S ) denotes the set of linear operators from TL® S 
to H® s . Let A G L(H® S ) be an arbitrary operator and S C {1, . . . , n} 
with | S"| = s. We denote by A[S] G L(H® n ) the embedding of the 
operator A into the Hilbert space 7i® n , i.e., the operator that acts as 
A on the qubits specified by S. 

An operator H : 'H® n — > n® n [ s called an s-local Hamiltonian if it 
is expressible in the form 

ff = I>iFi]. ( 2 ) 

j 

where each term Hj G Lffi®^^) is a Hermitian operator acting on a 
set Sj ■ | Sj | ^ s . 

A Hamiltonian is local if it can be expressed as a sum of terms, 
where each term acts on a bounded number of qubits. Indeed, in this 
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case, the corresponding time evolution can be approximately simu- 
lated by a universal quantum computer NCOO . 

For a direct physical implementation of the continuously varying 
Hamiltonian H(t) we require a stronger locality condition. Physical 
interactions are usually pair- interactions, unless one considers effective 
Hamiltonians. The system Hamiltonian can be thus decomposed as 

H = Y,H k i + Y, R k, (3) 

k<i k 

Hki is a Hermitian operator acting on the joint Hilbert space of particle 
k and I and is the free Hamiltonian of particle k. Furthermore, the 
interaction strength is decreasing with the distance. Therefore, we do 
not want to propose a scheme that relies on "weak" interaction terms 
among distant particles. We thus require that each particle is coupled 
to only a few other particles in its direct neighborhood. 

One of the most simple nontrivial examples are the Ising interac- 
tions on a 2D lattice. Our resource is the Ising Hamiltonian on an 
r x s rectangular lattices, i.e., 

Rising = £ 4 fc V«, (4) 
(k,l)eL 

where L are the edges of a rectangular lattice, i.e, a graph of order rs 
obtained by placing vertices at the coordinates | < i < r, < 

s < J} with edges joining just the pairs at unit distance. 

Let L' be a subgraph of L. We construct a final Hamiltonian 

Hp= Yl Wkl a z k)a z ] + Yl Wk(7 z k) with Wkl > w k^ Z > ( 5 ) 
(k,l)eL' k 

such that its ground states encode the solution to the NP-complete 
problem "independent set". Clearly, such Hamiltonians satisfy the 
locality condition. The aim of our paper is to show how such Hamil- 
tonians can be constructed using planar orthogonal embeddings of 
graphs and how they can be obtained efficiently from the 2D Ising 
model Hamiltonian -ffi s i ng . Together with the choice of a local initial 
Hamiltonian 

HB = Y4 k) (6) 

k 

our results allow to simulate efficiently the adiabatic quantum evolu- 
tion according to 

H(t) = {l-^)H B + ^H P . 
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2 Independent set problem 



The idea to consider the INDEPENDENT SET problem is motived by 
[EM EMI]. The INDEPENDENT SET problem |UT79] is defined 
as follows: 

• INSTANCE: Graph G = (V,E), positive integer v < \V\. 

• QUESTION: Does G contain an independent set whose cardi- 
nality is at least v, i.e., a subset V' C V such that \V'\ > v and 
such that no two vertices in V' are joined by an edge in El 

The INDEPENDENT SET problem remains NP-complete for cubic 
planar graphs GJS76 . A graph is called cubic if all vertices have 
degree 3, i.e., all vertices are connected to exactly three vertices. A 
graph is called planar if it can be drawn in the plane such that the 
edges do not intersect. An example of a planar cubic graph is shown 
in Figure 1. 



We consider a method that determines the maximal cardinality of 
independent set of a planar cubic graph. Let us recall how the solution 
to the maximum independent set can be encoded in the ground states 
of a pair-interaction Hamiltonian Hp. 

Theorem 1 (Planar spin glass within a magnetic field) 

Let G = (V, E) be a cubic planar. Determining the energy of the 
ground states of the corresponding Hamiltonian 
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Figure 1: Planar cubic graph 
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is equivalent to determining the maximum cardinality of independent 
sets of G. 

Proof. This has been shown in Bar82j (see also WB03 ). We include 
the proof here for completeness. We associate a variable Xk E {0, 1} 
to each vertex k € V. There is an independent set whose cardinality 
is at least v if and only if there is an assignment to the variables 
{Xk | k £ V} such that 

L = Y,X k - Y, Wi>v. (8) 

k£V (k,l)eE 

This is seen as follows. If V' is an independent set whose cardinality 
is at least v, then the assignment X k = 1 for k £ V' and Xk = for 
k £ V \ V fulfills inequality ©. 

Now let X±, . . . ,X n be an assignment that fulfills inequality (J5J). 
If V' = {k | Xk = 1} is not an independent set, then we must have 
\V'\ > v + p, where p := J2(k l)eE XkX\ > is the "penalty" for V 
not being an independent set. Let (k,l) € E with X^ = X^ = 1. By 
removing k from V 1 (i.e. setting X^ := 0) the cardinality of V drops 
by 1, while p drops by at least 1. After repeating this several times, 
we end up with an independent set whose cardinality is at least v. 

Setting Sk = 2X k - 1 for all k G V and observing that \E\ = ||V| 
for all cubic graphs, we obtain 

L = -\E S *-\ E S k S l + l\V\. (9) 
fcev (k,l)eE 

For E = —4L+\\V\ we see that there exists an independent set whose 
cardinality is at least k if and only if there is an assignment of values 
to the variables S k € {—1,1} (corresponding to the eigenvalues of a z ) 
such that 

E=J2$k+ E S k Si<^\V\-4v. (10) 

k£V (k,l)eE 

Now it is clear that determining the minimal energy E is equivalent 
to determining the maximal cardinality v of independent sets of G. □ 

In adiabatic quantum computing the initial Hamiltonian is chosen 

as 

h b = J2 ^ ( u ) 

fcev 
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and the time-dependent Hamiltonian as 

H{t) = (l-±)H B + ±Hp. (12) 

This Hamiltonian Hp does not necessarily satisfy the locality condi- 
tions. 

3 "Planar orthogonal" Hamiltonians 

Due to the lattice structure of our resource Hamiltonian we need to 
embed our graph into this structure. This can be done using planar 
orthogonal embeddings of graphs |KW01j . This idea is inspired by 
KL02 ( RAS02]. We shall be concerned with embedding graphs into a 
2D rectangular lattice. 

Definition 1 (Planar orthogonal embedding) 

A planar orthogonal embedding V of a graph G = (V, E) is a mapping 
that 

• maps vertices k to lattice points T(k) and 

• edges (k, I) 6 E to paths in the lattice such that the images of 
their endpoints T(k) and T(l) are connected and such that the 
paths do not share any vertices (besides the endpoints). 

Note that the map inserts "dummy vertices" if necessary to create 
the paths connecting the vertices and T[. A planar orthogonal 
embedding is shown in Figure 2. 

Every planar graph with maximum degree 3 admits a planar or- 
thogonal embedding on an \n/2\ x |_ n /2j- The algorithm presented 
in |Kan96j computes efficiently such planar orthogonal embeddings of 
graphs. We used AGD (Libary of Algorithms for Graph Drawing) to 
compute the embedding |AG D02 . 

In the proposal of KL02 the Hamiltonian Hp is considered. The 
planar orthogonal embedding gives a regular wiring among the qubits. 
This means that the couplings are not spatially local. In contrast, we 
need a Hamiltonian Hp that contains only nearest-neighbor interac- 
tions. This is necessary that it can be simulated by -fusing- The 
idea is to use the dummy vertices as wires that propagate the state 
of a (real) vertex spin to the neighborhood of another vertex. This 
can be achieved by constructing a path of adjacent dummy vertices, 
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Figure 2: Planar orthogonal embedding of the graph in Fig. 1 

each interacting with its neighbor by a strong ferromagnetic coupling. 
Furthermore, the first dummy at one end of this "dummy path" is 
strongly ferromagnetically coupled to a vertex and the last dummy 
at the other end is in the neighborhood of another real vertex, cou- 
pled to it via a usual antiferromagnetic interaction. The interaction 
strength is chosen in such a way that it is always energetically better 
when all dummies have the same state as the real vertex to which they 
are connected to than to have a mismatch along the "ferromagnetic 
path" . 

Formally, this construction is as follows: 

• The dummy vertices have no local a z term. 

• The vertices T(k) have a z as local Hamiltonians. 

• Let (k, I) € E be an edge of G. 

If Tfc and T; are adjacent, then the coupling between and T[ 
is chosen to be antiferromagnetic, i.e., a z ® a z . 

Otherwise there are m dummy vertices v% , . . . , v m such that the 
path (Tk,vx, . . . ,v m ,Ti) connects the vertices and T;. The 
couplings between and v\ and Vi and Vi + \ for i = 1, . . . , m — 1 
are chosen to be ferromagnetic with coupling strength c, i.e., 
— ca z (8> cr z . The coupling between v m and Ti is chosen to be 
antiferromagnetic, i.e., a z <g) a z . 
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The corresponding "planar orthogonal Hamiltonian" is shown in Fig- 
ure 3. The filled circles correspond to dummy vertices that do not 
have any local Hamiltonian. The circles with indices correspond to the 
original vertices of G. They have a z as local Hamiltonians. The thin 
lines correspond antiferromagnetic interactions and the thick lines to 
ferromagnetic interactions. The idea behind this construction is that 




Figure 3: Hamiltonian corresponding to the planar orthogonal em- 
bedding in Fig. 2 

there is a direct one-to-one correspondence between the ground states 
of Hp and Hp. The same is true for the first excited states. This can 
be seen as follows: 

Let (k, I) G E be an edge of G and v±, . . . , v m ,Ti) be the path 
on the lattice connecting and T/. The variables Sr k , Si, ... , S m £ 
{0, 1} indicate whether the corresponding qubit is spin up or spin 
down. 

A ground state satisfies the condition that S\ , . . . , S m are all equal 
to Sr k - To show this we define the number of mismatches along the 
path to be the number of occurrences of 5r fe ^ Si, Sj ^ Sj + i for 
i = 1, . . . ,m — 1. This number is denoted by 5. 

Then the minimal possible energy (due to the couplings along the 
path) is 

c{-m + S)-l. (13) 



<s 



If we remove the mismatches (by setting S% := SW for i = 1, . . . , m) 
then the maximal possible energy is 



cm + 1 . 
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By choosing c = 3 minimal energy can be achieved only if the states of 
all dummy vertices are equal to the state of the qubit corresponding 
to T k . 

For adiabatic quantum computing it is important that the gap 
between the ground and first excited states of the Hamiltonian at all 
times is sufficiently large. We show that the modification of Hp to 
Hp does not decrease this gap. 

The gap between the ground and the first excited states of Hp 
is smaller or equal to 8. This is seen as follows. Let Si,...,S n G 
{— 1,+1} be an assignment corresponding to a ground state of Hp. 
Pick any vertex k and let l%,l2,h be at the three vertices connected 
to k. By flipping S*. the energy can increase by at most 8 because the 
relevant Hamiltonian is 



By choosing c = 9 it is seen that the first excited states of Hp satisfy 
the condition that the states all of dummy vertices are equal to the 
vertex of IV 

4 Simulating "planar orthogonal" 
Hamiltonians 

To implement the time-evolution according to the Hamiltonian Hp 
we make use of the concepts of simulating Hamiltonians that has been 
used in nuclear magnetic resonance for a long time EBW87 . These 
techniques rely on the so-called average Hamiltonian approach. The 
idea is to conjugate the natural time evolution by unitary control 
operations Uj, i.e., the total evolution is 

ttfc exp(— iHtk)v/ k . . .u\ ex.Y>{—iHt2)u\ui exp(— iHt\)u[ , 

where the system evolves in an undisturbed way during periods of 
length ti,t2, ■ ■ ■ ,tk- If these periods are short compared to the time 



3 



a 




i=i 



9 



scale of the natural evolution, the total dynamics is approximatively 
the same as the evolution according to the average Hamiltonian 

j 

with t := Usually, the control operations on n particles are 

assumed to be of the form 

u := v\ <g) V2 <8> • • • <S> v n 

where Vj is a unitary acting on particle j. The design of simula- 
tion schemes for Hamiltonians with n particles interacting via pair- 
interactions leads to non-trivial combinatorial problems (e.g. |LCYY00l 
IDlNBT()lllWJB()2[ITWr302] b An experiment al proposa l for simulating 
dynamics in optical systems is presented in |JVD + 02] . 

Starting from the Ising Hamiltonian Rising, we can implement the 
Hamiltonian Hp with time overhead (slow-down) 2c + 1 and 16 time 
steps by interspersing the time evolution according to -fusing by local 
operations in X (g> X <g) • • • (g> X, where X = {1, a x }. 

Following the ideas of |LCYY0f)l IW.7B02| we construct a selective 
decoupling scheme based on Hadamard matrices. Due to the special 
form of -ffising it is sufficient to use the Hadamard matrix of size 4 
only. 

Our scheme consists of 4 subroutines that implement the following 
couplings of Hp: 

1. horizontal a z <X> a z , 

2. vertical a z ® a z , 

3. horizontal —ca z <8> cr z , and 

4. vertical — ca z (g> a z 

The indices i,j enumerate the rows and the columns of the lattice. 
We denote the columns of the Hadamard matrix of size 4 

1 1 1 \ 

1 1 -1 

1 -1 -1 
1-1 I ) 

by W{0, 0), W(0, 1), W(l, 0) and W(l, 1). 

Let v = (v\,V2,V3,V4 : ) € {—1, l} 4 be a column vector. We use the 
abbreviation 



W :-- 
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"apply v at 

to denote the following control sequence with 4 equally long time steps: 
at the beginning and the end of the sth time step we apply a x on the 
qubit at position if v s = —1 and do nothing if v s = 1, where time 
step s runs from 1 to 4. 

Let v,v' G {— 1,1} 4 . One easily checks that applying v and v' at 
adjacent lattice points changes a z <8> a z to (v, v') a z (g> a z , where (v, v') 
denotes the inner product of v and v'. This is the key observation for 
constructing the selective decoupling scheme. 

In the first and second subroutines the length of the 4 time steps is 
chosen to 1/4. Let us consider the first subroutine. The vertical cou- 
plings are automatically removed if we apply in rows with even indices 
only W(0, 0) and W(1,0) and in rows with odd indices W(1,0) and 
W(l, 1). The choice between W(a, 0) and W(a, 1) depends on which 
horizontal interactions should remain or be switched off. Explicitly, 
this choice is as follows. Choose W(a,0) for the leftmost spin. If the 
interaction between the spins (j — 1) and j should remain, then apply 
the same W(a, b) to j as to (j — 1). Otherwise (i.e. the coupling should 
be switched off) apply the second possible W(a, b') to j. 

The second subroutine is obtained from the first subroutine by 
exchanging the roles of rows and columns of the lattice. 

In the third and fourth subroutines the length of the 4 time steps 
is chosen to c/4. The third subroutine is obtained from the first sub- 
routine by apply (— l) 3 v instead of v to the spin j. Finally, the fourth 
subroutine is obtained from the third subroutine by exchanging the 
roles rows and columns of the lattice. 
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